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Abstract 

We generalize the virial theorem to the warped DGP brane-world scenario and consider its im- 
plications on the virail mass. In this theory the four-dimensional scalar curvature term is included 
in the bulk action and the resulting four-dimensional effective Einstein equation is augmented 
with extra terms which can be interpreted as geometrical mass, contributing to the gravitational 
energy. Estimating the geometrical mass M{r) using the observational data, we show that these 
geometric terms may account for the virial mass discrepancy in clusters of galaxies. Finally, we 
obtain the radial velocity dispersion of galaxy clusters ay (r) and show that it is compatible with 
the radial velocity dispersion profile of such clusters. 

PACS numbers: 04.50.-h, 04.20.Jb, 04.20.Cv, 95.35. +d 
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1 Introduction 

In recent times theories with higher-dimensions have become popular in high energy physics, specially 
in the context of hierarchy problem and cosmology [1]. In this scenario it is purported that our four- 
dimensional universe is a subspace called the brane, embedded in a higher-dimensional space-time 
called the bulk. One of the most successful of such higher-dimensional models is that proposed by 
Randall and Sundrum whose the bulk has the geometry of an AdS space admitting Z2 symmetry [2] . 
They were successful in explaining what is known as the hierarchy problem: the enormous disparity 
between the strength of the fundamental forces. The Randall Sundrum (RS) scenario has had a great 
impact on our undersetting of the universe and has brought higher-dimensional gravitational theories 
to the fore. In certain RS type models, all matter and gauge interactions reside on the brane while 
gravity can propagate into the bulk. Using the Israel junction conditions [3] and the Gauss-Codazzi 
equations, one can obtain the field equations on the brane, as employed by Shiromizu, Maeda and 
Sasaki (SMS) [4]. There are two very important results that arise from the effective four-dimensional 
Einstein equations on the brane. The first one is quadratic energy-momentum tensor, tt^, which is 
relevant in high energy and the second one is the projected Weyl tensor, £^ u , on the brane which is 
responsible for carrying on the brane the contribution of the bulk gravitational field. The cosmological 
evolution of such a brane universe has been extensively investigated and effects such as a quadratic 
density term in the Friedmann equations have been found [5, 6, 7]. 

An alternative scenario was subsequently proposed by Dvali, Gabadadze and Porrati (DGP) [8]. 
The DGP proposal rests on the key assumption of the presence of a four-dimensional Ricci scalar in the 
bulk action. There are two main reasons that make this model phenomenologically appealing. First, 
it predicts that four-dimensional Newtonian gravity on a brane-world is regained at distances shorter 
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than a given crossover scale r c (high energy limit), whereas five-dimensional effects become manifest 
above that scale (low energy limit) [9]. Second, the model can explain late-time acceleration without 
having to invoke a cosmological constant or quintessential matter [10]. An extension of the DGP 
brane-world scenario have been constructed by Maeda, Mizuno and Torii which is the combination of 
the RS II model and DGP model [11]. In this combination, an induced curvature term appears on the 
brane in the RS II model. This model has been called the warped DGP brane-world in literature [12]. 
In this paper, we consider the effective gravitational field equations within the context of the warped 
DGP brane-world model and obtain the spherically symmetric equations in this scenario. So much 
for the success of the DGP model, a word of caution is in order; the theory predicts the existence 
of ghost-like excitations. Many scenarios have been undertaken to explain away such ghosts, but as 
yet no satisfactory solution exists. The interested reader should consult [13] for further insight. We 
do not discuss such excitations since our aim lies in studying the virial mass discrepancy in warped 
DGP models. 

Modern astrophysical and cosmological models are faced with two severe theoretical difficulties 
which can be summarized as dark energy and dark matter problems. The problem of dark matter 
is a longstanding problem in modern astrophysics. Two important observational issues, the behavior 
of the galactic rotation curves and the mass discrepancy in clusters of galaxies led to the necessity 
of considering the existence of dark matter at the galactic and extra-galactic scales [14]. The total 
mass of a cluster can be estimated in two ways. One can apply the virial theorem to estimate the 
total dynamic mass My of a rich galaxy cluster from measurements of the velocities of the member 
galaxies and the cluster radius from the volume they occupy. The second is obtained by separately 
estimating the mass of each individual members and summing them up to give a total baryonic 
mass M. It is always found that My is greater than M. This is known as the missing mass 
problem. Nevertheless, the existence of dark matter was not firmly established until time when the 
measurement of the rotational velocity of stars and gas orbiting at a distance r from the galactic 
center was performed. Observations show that the rotational velocity increase near the center of 
galaxy and remain nearly constant. This discrepancy between the observed rotation velocity curves 
and the theoretical prediction from Newtonian mechanics is known as galactic rotation curves 
problem (Figure 1). These discrepancies are explained by postulating that every galaxy and cluster 
of galaxy is embedded in a halo made up of some dark matter [14]. To deal with the question of 
dark matter, a great number of efforts has been concentrated on various modifications to the Einstein 
and the Newtonian gravity [15, 16, 17, 18, 19]. Several theoretical models, based on a modification 
of Newton's law or of general relativity, have been proposed to explain the behavior of the galactic 
rotation curves. In the modified Newtonian dynamics (MOND) theory which has been proposed by 
Milgrom [20], the Poisson equation for the gravitational potential, V 2 </> = AirGp, is replaced by an 
equation of the form V[/i(x)(|V^|/ao)] = 4irGp, where ao is a fixed constant and a{x) a function 
satisfying the conditions u(x) = x for x « 1 and u(x) = 1 for x » 1. The force law, giving the 
acceleration a of a test particle becomes a = ajy for » ao an d a = y / a/yao f° r &N « ao, where 
cijv is the usual Newtonian acceleration. The rotation curves of the galaxies are predicted to be flat, 
and they can be calculated once the distribution of the baryonic matter is known. A relativistic 
MOND inspired theory was developed by Bekenstein [21]. In this theory gravitation is mediated 
by a metric, a scalar field and a four-vector field, all three dynamical. For alternative theoretical 
models to explain the galactic rotation curves, see [22]. One other such modification is that of the 
RS brane-world scenario [23]. It has been argued that a modified theory of gravity based on the 
RS brane-world scenario can explain the observations of the galactic rotation curve of spiral galaxies 
and the virial theorem mass discrepancy in clusters of galaxies without introducing any additional 
hypothesis [23]. 

Our main purpose in this paper is to obtain the generalized form of the virial theorem in Warped 
DGP brane-world model by using the collision-less Boltzmann equation. In what follows, we first give 
a brief review of warped DGP brane-world model and the gravitational field equations are derived 
in this model. Next, we use the relativistic Boltzmann equation to derive the virial theorem which 
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Figure 1: Observed rotation velocity curve of the NGC3198 (dotted line) and the prediction from Newtonian theory 
(solid line). 



is modified by an extra term which may be used to explain the virial mass discrepancy in clusters 
of galaxies. Finally, we identify the geometrical mass of cluster in terms of the observable quantities 
and obtain the radial velocity dispersion of galaxy clusters. 

2 Effective field equations on warped DGP brane 

Let us start by presenting the model used in our calculation [11]. Consider a five-dimensional space- 
time with a four-dimensional brane is located at Y(X A ) = 0, where X A , (A = 0,1,2,3,4) are 
five-dimensional coordinates. The effective action is given by 



S — Sbulk + Sbra 



where 



•Sbulk 



I 



24 



(2) 



and 



Sbr 



J 

JY=0 



(3) 



where 



8-irGx, is the five-dimensional gravitational constant, 1Z and C m are the five-dimensional 
scalar curvature and the matter Lagrangian in the bulk, respectively. Also, x^, (/j, = 0, 1, 2, 3) are the 
induced four-dimensional coordinates on the brane, is the trace of extrinsic curvature on either 
side of the brane [24, 25] and A>mne(<7a/35 ^) is the effective four-dimensional Lagrangian, which is 
given by a generic functional of the brane metric g a p and matter fields. 
The five-dimensional Einstein field equations are given by 



k 2 



Tf B + 5(Y) TAB 



(4) 



where 



T (5) _ 
1 AB = 



5C 



(5) 



'sg AB 



(5) 



and 



ry ^^brane ^ 

= ~^ ~ : r g^iuJ-brane- 



5g> 



(6) 



We study the case with an induced gravity on the brane due to quantum corrections [8]. The 
interaction between bulk gravity and the matter on the brane induces gravity on the brane through 
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its quantum effects. If we take into account quantum effects of matter fields confined on the brane, 
the gravitational action on the brane is modified as 



£brane(9a/3, VO — ^R ~ + An, 



(7) 



where fi is a mass scale which may correspond to the four-dimensional Planck mass, A& is the tension 
of the brane and An presents the Lagrangian of the matter fields on the brane. We note that for 
Xb = and A^ 5 ) = action (1) gives the DGP model and gives the RS II model if fi = 0. 
We obtain the gravitational field equations on the brane-world as [4] 



Gfiu — 



2k{ 



(8) 



V„t£ = -2Tf B n A g*, (9) 
where V v is the covariant derivative with respect to and the quadratic correction has the form 

(10) 



7T, 



1 



24 



and the projection of the bulk Weyl tensor to the surface orthogonal to n A is given by 

p _ r (5) A B C D 



fill 



The symmetry properties of imply that in general we can decompose it irreducibly with respect 
to a chosen 4-velocity field v M as [6] 



£\xv = — I — I 



U (v^Vv + -h^) + 2Q (p v u) + 



(12) 



where h^ u = g^ u + v^v u projects orthogonal to and the factor (K5//C4) is introduced for dimension 
reasons. Here 

is an effective non-local energy density or "dark radiation" term on the brane, arising from the free 
gravitational field in the bulk, = (K4/K5) 4 h°£ a ^v^ is an effective non-local energy flux and 



U = - 



111/ 



h t . a h 



(P v) 



£ 



a/3, 



is a spatial, symmetric and trace-free tensor. In what follows for the static spherically symmetric 
brane we have Q u = 0, and we may choose P^ v = P{r){r^r u — 5/1^), where the "dark pressure" P(r) 
is a scalar functions of the radial distance r, r M is a unit radial vector and at any point on the brane 
in inertial frame = 5q, h^ v = diag(0, 1, 1, 1) [26]. 

In order to find the basic field equations on the brane with induced gravity, we have to obtain the 
energy- momentum tensor of the brane t^ u , given by definition (6) from the Lagrangian (7), yielding 



(13) 



Assuming that the five-dimensional bulk space includes only a cosmological constant A^ 5 ^ and inserting 
equation (13) into equation (8), we find the effective field equations for four-dimensional metric g^ v 
as 



1 + ^4fi 2 ) GV = ±\ b 4 V - A 4y > - 4(i 2 K^ a pG a P + 4 Trg? + A£? 



1 



(14) 
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where 



Kfiupa — . {dpL'^pu QppTpa Qva^pp 

71 p,u = ~^T^ a T u + —TT^v + —gp U T a pT @ — —gp V T , 



24^ 



7T 



(G) 
pu 



-^GpaG® + —GG^ + -gpvG al 3G al3 - -^gpvG 2 , 



and the effective cosmological constant on the brane is given by 



1 



A (5) + W,Xl 



6 -5-6 



(15) 
(16) 
(17) 

(18) 



We note that for p, = these equations are exactly the same effective equations as in the reference 
[4]- 



3 Field equations for a cluster including identical and collision-less 
point particles 

Now we consider an isolated and spherically symmetric cluster being described by a static and spher- 
ically symmetric metric 

ds 2 = -e x(r) dt 2 + e u{r) dr 2 + r 2 {d9 2 + sin 2 9d<p 2 ). (19) 

Suppose that the clusters are constructed from identical and collision-less point particles (galaxies). 
This multi-particle system can be described by a continuous non-negative function /s(x M ,p M ), dis- 
tribution function, which is defined over the phase space. In terms of the distribution function the 
energy-momentum tensor can be written as [27, 28] 



Tp U = J fBmVfj.Vudv, (20) 

where m is the mass of each galaxy, is the four-velocity of the galaxy and dv = ^dvrdvgdv^, is the 
invariant volume element of the velocity space. We also assume that the matter content of the bulk 
is just a cosmological constant A^ 5 ) and the energy-momentum tensor of the matter in a cluster of 
galaxies can be represented in terms of spherically symmetric perfect fluid as 

Tpv = {Pb + Pb) VpV v + n9pv, (21) 
where v^v^ = —1. Use of equations (20) and (21) leads to the following relation for pb and pb 

Pb = p(v 2 ), Pb = p(v 2 ) = p{vj) = p(vl) , (22) 

here (v 2 ) represents the usual macroscopic averaging which is defined as (v 2 ) = - J fBV 2 mdv where 
p is the mass density and in term of Jb is given by p = J fsmdv [29, 30]. 

Using equation (14), the gravitational field equations on the brane become 

(1 + ^4" 2 )^ ("I + rv + J) = \\b4pb + A 4 + 4p 2 n + %l - <P A A {G) + 4)~U(r), (23) 



12'° ° r u K\\ h 

1 + y4" 2 )^(l + rA'-e*) = iA 6 ^ 
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(1 + ^ K iii 2 ) e — (2A' - 2i/' - X'u'r + 2A"r + X' 2 r) = ^X b 4p b - 2A 4 - 2^ 2 K 2 + ^ (p 2 . + 2p 6 p 6 ) 



K5 
6 



+ 24^f G) + ^-[U(r)-P(r)l 



K 2 A X b ' 



(25) 



where U(r) and P (r) are the dark radiation and dark pressure and K^^G^ = K.^ v term is given 
by 

\ {Pb + Pb) + VbX'v' - 2p b X" - p b X' 2 + \{p b + p b )e u 



{pb + Pb) ~ k>b 



r 



e v (p b - 3p b ) 



24 



2e u 2u' 2 AX' w , w . w2 
+ ? + AV - 2A" - A 

ry£ ry ryZi ry 



+ 



Pbe 



24 



4z/ 4A' 



— - — - 4 + AV - 2A" - A' 2 + ^~ 



4e^ 



(26) 



= -V 



2z/ 



pAV - 2pA" - P X' 2 { Pb - p b ) + ^(p b + p b ) e v + -r ( Pb + p b ) 



e u (p b - 3p b ) 



24 



4u' 2e v 



+ 



2A' 



^2 ^2 



+ AV - 2A" - A' 2 



Pbe' 



24 



4z/ 46^ 4 

+ 



4A' 



+ AV - 2A" - A' 2 



ry ^*2 ^*2 ^* 



(27) 



r 2 _ r3 _ 
^2 — ^3 — 



PbX' pb^_ {pb ~ Pb) e u _ (p b - Pb) 



e v (p b - 3pb) 



48 



61/ 8 6A' w , nW/ x/2 
— + ^ ^ + AV - 2A" - A' 2 



ry ty^t ty^i ty 



Pbe 



24 



4z/ 4e 



4A' 



— + ^—--- — + AV - 2A" - A' 2 



and use of equation (17) leads to the components of tt^ G ^ as 



7T, 



0(G) 

— 



7T 



1(G) 



-6z/ 



24 



2A AV 2A A 2 \ e 
+ -5- + — I + 



-2v 



T 
1 „ // 



2 /f->2 / 192 ^3 



+ 



4A A v 
16A 



8A A 



+ 



4AV 



+ A 2 i/ 2 -2AV +4A"A' : 



+ 4A" 2 + A' 4 - % 



-r v ( 2v 



16A' 


81/' A 




r 


\2v' 2 


4A' 3 




r 



(28) 



+ 



16A 1/ 



(29) 



24 



^3 ^>3 ^->2 ^2 ^2 / ]_92 



6X_AV 2X_ A^\ ; r 



2 ^ 16z/_48V 8z/ A" 16AV 



4A 1/ A + 



8A'A" 4AV 2 



+ //V 2 -2A'V +4A'A' 2 + 



16A //n 

— ^ + 4A 2 + A 4 - 



r r 

20A' 2 , 



4z/ 2 4A' 3 



+ 



(30) 



7T. 



2(G) 



7T 



3(G) 



-( 


4 AV 


_ 2A" A' 2 


24 \ 


/y2 ry2 


2A'A" 


AV 2 


2v 2 A' 3 



2e v 



-2u 



48 



2A 2 i/ 2A 1/ 



4A" 



+ ^— + — Z1 2- + — + ^r), (3i) 

ry£t fy fy ry £t ry fyii 

where a prime represents differentiation with respect to r. In the next section, we will investigate the 
influence of the bulk effects on the dynamics of the galaxies in warped DGP brane-world model. 
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4 The virial theorem in warped DGP brane 

In order to derive the virial theorem for galaxy clusters, we have to first write down the general 
relativistic Boltzmann equation governing the evolution of the distribution function fg. The galaxies, 
which are treated as identical and collision-less point particles, are described by this distribution 
function. For the static spherically symmetric metric given by equation (19) we introduce the following 
frame of orthonormal vectors [27, 28, 29] 



(i) 
p 



-7> 



er y = r8 , 



&Z' = t sin 9Sg, 



where g^effiefi 



(a)(6) _ -p^g f our _ velocity of a typical galaxy with 



vv,, 



components is written as 



,(«) 



Ve 



a = 0,1,2,3. 



The relativistic Boltzmann equation in tetrad components is given by 

Xa)o 9f B 



(a) Q x p '(6)(c) 



(32) 

— 1, in tetrad 
(33) 

(34) 



where fs = fB{x^,v^) and 7(&v c ) = e n& e (b) e (c) are ^ ne distribution function and the Ricci rota- 
tion coefficients, respectively. Assuming that the distribution function is only a function of r, the 
relativistic Boltzmann equation becomes 



dfs 
dr 



v 2 8\ {v 2 6 + vp 
2 dr r 



dfs _ Vr_ 

dv r r 



df B ^ df B 



dv u 



W2, 



cote 



dfB df B 



dv u 



dvg 



0,(35) 



where we have defined 



v t , v {1) = v r , 



■U (2) = Vg, 



(36) 



Since we have assumed the system to be spherically symmetric, the term proportional to cot# must be 
zero. Multiplying equation (35) by mv r dv where dv = ^dvrdvgdv^, and integrating over the velocity 
space and assuming that the distribution function vanishes rapidly as the velocities tend to ±oo, we 
obtain 



d_ 

dr 



P(v 2 r) 



1 

+ 2 P 



<«?> + {v 2 ) 



dX 
r dr-~ P 



(v 2 ) + {vp - 2{v 2 r ) 



0, 



(37) 



where p is the mass density and (v 2 ) represents the average value of v 2 . Multiplying equation (37) 



by 47rr and integrating over the cluster leads to 

rR 



j o '47rp[(v 2 r ) + (v 2 e ) + (vP 
This equation is reduced to 



r 2 dr 



1 rR 
~ 2 Jo 



47rr 3 p (v 2 ) + {v 2 ) 



dX , 

— dr = 0. 

or 



2K 



\j\^p[{v 2 ) + {v 2 .)\ ^dr = 0, 



2 Jo -■ r ' N r/ J dr 

where the total kinetic energy of the galaxies is defined as 

rR 



K = J q 2irp [(v 2 ) + (v 2 e ) + (v 2 p 



r 2 dr. 



(38) 



(39) 



(40) 



Now, using these relations and adding the gravitational field equations (23)- (25) and equation (22) 
we find 
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1 i ^ b 4 2 , 



r 



aV a" a' 2 \ «l 

— + t + - = f p 



<«?> + + (vl) + (vi) 



+ 



2\ b F 

4 4 



<^ 2 > 2 + <*v 2 > 2 + (v 2 ) 2 + (4)'< 



0(G) 



K4A5 



4 2 

- A4 H — [ICq — JC{- 2/C 2 ] 



I7(r). 



(41) 



where |AfeK| = In order to obtain the generalized virial theorem we have to use some ap- 
proximations. First, since the dispersion of the velocity of galaxies in the clusters is of the order 



600 - WOOkm/s, i.e., ( 



v\2 



4x10 — 1.11 x 10 < 1, we can neglect the relativistic effects in 



the relativistic Boltzmann equation and use the small velocity limit approximation. In other words, 
(v 2 ) ps (vq) ps (t> 2 ) <C (v 2 ) ss 1. Second, the intensity of the gravitational effects can be estimated 
from the ratio GM/R, which for typical clusters is of the order of 10 -6 <C 1. Therefore inside the 
galactic clusters the gravitational field is weak and we can use the weak gravitational field approxi- 
mation. Then the term proportional to AV and A' 2 in equation (41) may be ignored. Thus, assuming 
that e A ps e u ps 1 inside the cluster [23], we can write equations (41) as 



(1 + 



A;, 4 2 



K4 



tf) + (V 2 r ) + (vj) + (Vl) 



A 4 



K 4 I, 2 

6r 



2,v'(v 2 )+2\'(v 2 )+\"r(v 2 ) + Zu'(v 2 r ) 



+ 



6 



k 4 A& 



I7(r). 



(42) 



On the other hand, for clusters of galaxies the ratio of the matter density and of the brane tension is 
much smaller than 1, p/\ << 1, so that one can neglect the quadratic term in the matter density in 
above equation. These conditions certainly apply to test particles in stable circular motion around 
galaxies, and to the galactic clusters. Thus, we can rewrite equation (42) as 



! , X b 4 2\ 1 9 
1 + ^ 



where 



dr 



K 2 



p-A4 + 4fi 2 [V(r) + 3U(r)] + 



6 



K4A5 



U(r) 



U(r) 



.EL 

6r ' 



■p( r \ = — JL r 2 — 

6r 2 dr \ dr 

Multiplying equation (43) by r 2 and integrating from to r yields 



(43) 



K 2 



1 



— M(r) + -A 4 r 3 



-M DGP {r) 



K4. 



tor'"*"" v ' 8vr 

The total baryonic mass and the geometrical masses of the system are given by 



M RS (r) = 0. 



M(r) = 4tt 



I p( r ') r 2 d r > 
Jo 



and 



4Mrs(t) = ^- fu(r')r' 2 dr\ 
K4A6 jo 



(44) 



(45) 



(46) 



4M D Gp(r) = 8tt4h 2 [ r [P(r')+3U(r')]r' 2 dr', 
Jo 



(47) 
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where 

M{r) = Mrs (r) + M DG p (r) . (48) 

Multiplying equation (44) by dM ^ an d integrating from to R, we finally obtain the generalized 
virial theorem in warped DGP scenario as 

4 2^7/ TXT" , 1 

where 



1 + y 4n )2K + W + -A 4 I + W = 0, (49) 



W 
and 



W = -^i f R Mid d M(r), (50) 

87T Jo f 
K 2 /■• R 

— - / M(r)prdr, (51) 
2 Jo 



R 

r 2 dM(r), (52) 



where is the gravitational potential energy of the system. At this point it is worth nothing that 
for p = 0, we have Mdgp = and the virial theorem in the warped DGP brane- world is reduced to 
the virial theorem in the RS brane scenario [23] 

2K + W + \kj + Wrs = 0, (53) 



where 



k 2 rR 



W RS = -^ M RS (r)prdr. (54) 
2 Jo 

As one can see the gravitational energy modified by Wrs which has its origin in the global bulk effect 
due to the term. We can also recover the virial theorm in the standard general relativity from 
equation (53) in the limit A^ 1 — > 0. A alternative possibility in recovering the four-dimensional virial 
theorem is to take the limit K5 — > 0, while keeping the Newtonian gravitational constant n\ finite 
[11]. 

In the case A;, = A^ 5 ) = 0, the virial theorem in warped DGP brane-world is reduced to the virial 
theorem in DGP brane scenario [19] 

2K + W + W DGP = 0, (55) 

where 

k 2 f R 

W DGP = -^ M DGP (r)prdr. (56) 
2 Jo 

Note for a Minkowski DGP bulk space we have = 0, thus in above equation Wrs = 0. There 
is difference between our model with references [23] and [19]. The virial theorem in warped DGP 
brane-world modified by both Wrs and Wdgp-, which the first is due to the global bulk effect whereas 
the second term has its origins in the induced gravity on the brane due to quantum correction. 
Now, we introduce the radii Ry , Ri and TZ as 

M 2 (r) 
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o f*r 2 dM(r) , s 



K- KlM2{r) (59) 



where Ry is the virial radius and TZ is defined as the geometrical radius of the clusters of galaxies. 
Defining the virial mass as [29] 



2K = -^^-, (60 ) 



8tt R v 

and using the following relations 



- 2 M 2 



w = -^in I = MRl (61) 



8vr R v 

the generalized virial theorem (49) is simplified as 

(, M A b 4 2\ (M V \ 2 1 8nA 4 R V R 2 ( M\ 2 ( Ry\ 

We have three types of mass in equation (62), namely, the total baryonic mass of the system repre- 
sented by M (including the baryonic mass of the intra-cluster gas and of the stars, other particles like 
massive neutrinos), the virial mass represented by My and finally, the geometrical mass represented 
by M. 

On large distance scales associated with galaxies, we can ignore the contribution of the effective 
cosmological constant to the mass energy of the galaxy. Also, it is found that My is considerably 
greater than M for most of the clusters and we can neglect the term unitary in equation (62). 
Therefore, the virial mass is given by 

My{r)^M{r)[^ 1/2 . (63) 
This equation shows that the virial mass is proportional to the geometrical mass. 



5 Astrophysical applications 

In this section, we emphasize that the astrophysical observations together with the cosmological 
simulations have shown that the virilized part of cluster is a measure of a fixed density such as a 
critical density, p c (z) at a special red-shift, so that py = 3My /4-irRy = Sp c (z), where My and Ry are 
the virial mass and radius, respectively. As is well known, p c (z) = 1i 2 (z)3Hq /8ttG, where the Hubble 
parameter is normalized to its local value, i.e., h 2 (z) = f2 m (l + z) 3 + Oa, where Q m and S7a are the 
mass density and dark energy density parameters, respectively [31]. By knowing the integrated mass 
of the galaxy cluster as a function of the radius, one can estimate the appropriate physical radius for 
the mass measurement. The radii commonly used are r2oo or r 500- These radii lie within the radii 
corresponding to the mean gravitational mass density of the matter ptot = 200p c or 500p c . A useful 
radius is r2oo to find the virial mass. The numerical values of the radius r2oo for the cluster NGC 
4636 are in the ranges r2oo = 0.85 Mpc and r2oo = 4.49 Mpc for the cluster A2163, so one can deduce 
that a typical value for r2oo is 2 Mpc. The masses corresponding to r2oo and rsoo are denoted by 
M200 and M500, respectively, and it is usually assumed that My = M200 and Ry = T2qq [32]. 
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5.1 Geometrical mass estimated using the Jean's relation 

Now, we are going to obtain Ai(r) as a function of r by comparing the virial theorem results with 
the observational data for galaxy cluster which can be obtained from the X-ray observation of the 
gas in the cluster. The most of the baryonic mass in clusters is in the gas form, therefore we assume 
that the energy-density and pressure in T^ u is that of a gas as 

P = Pg(r), P = Pg(r)- (64) 

In majority of clusters most of the baryonic mass is in the form of the intra-cluster gas. The gas mass 
density p g {r) distribution can be fitted with the observational data by using the following expression 
for the radial baryonic mass distribution [32] 

-3/3 

Pg(r)=Po(l + ^j 2 , (65) 

where r c is the core radius, and po an d (3 are cluster independent constants. A static spherically 
symmetric system of collision-less particles that is in equilibrium, can be described by the Jean's 
equation [14] 

d r / n 2l , 2 Pg( r ) ( 2 2 \ / dMj) 

[Pg(r)<T r \ + — jj-— [a r - o^) = ~Pg( r )-^, ( 66 ) 

where $(r) is the gravitational potential, ay and <r| are the mass-weighted velocity dispersions in 
the radial and tangential directions. We assume that the gas is isotropically distributed inside the 
cluster, so a r = ag^. The gas pressure is related to the velocity dispersion and gas profile density by 
Pg = Pg a r • By assumption that the gravitational field is weak so that it satisfies the usual Poisson 
equation 2V 2< 1 ) ~ n\ptoti where ptot is the energy density including p g and other forms of matter, like 
luminous matter and the geometrical matter, etc., the Jean's equation becomes 

dp g (r) , ,dMr) 4 M tot , , f(t -. 

where Mt t(f) is the total mass inside the radius r. The observed X-ray emission from the hot ionized 
intra-cluster gas is usually interpreted by assuming that the gas is in isothermal equilibrium. Therefore 
we assume that the gas is in equilibrium state having the equation of state p g {r) = ^0^-pg{r), where 
kg is Boltzmann constant, T g is the gas temperature, p = 0.61 is the mean atomic weight of the 
particles in the gas cluster and m p is the proton mass. Equation (67) then gives 



pm p n\ dr 



M tot (r) = - -z^dtr 2 X ln P9 W ■ (68) 



Now, use of the density profile of the gas given by equation (65) leads to the mass profile inside the 
cluster as 

2^k B T g f3 r 3 

M tot {r) = 1 q—^- 69 

On the other hand, using equation (45) and (48) we can obtain another expression for the total mass 

dM tot (r) 2 48vr m * 2 , Sttk^ 2 r _, . . „ m 2 /-^ 
= 4-irr A p g (r) + -^—U(r)r^ H ^ — [P{r) + 3U(r)\ r z , (70) 

(XT ^4^6 ^4 

substituting equations (69) and (65) into equation (70) we obtain the following expression 

+ 2 4f + «M] - - ,0 f 1 + . (71) 



n\\ h k\ 1 ' pmpKl (r 2 + r 2 ) 2 
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Finally, substituting above equation into equation (48), in the limit r » r c considered here, we 
obtain the following geometrical mass 



M(r) ~ 



lirripKi c 3(1 - /3) 



(72) 



which includes both the local and non-local bulk effects. Observations show that the intra cluster 
gas has a small contribution to the total mass [31, 32, 33, 34], thus we can neglect the contribution 
of the gas to the geometrical mass and rewrite equations (72) as 



24^ \ 
^ fim p K 2 4 J 



Now Let us estimate the value of M.{r). First we note that ksTg « 5KeV for most clusters. The 
virial radius of the clusters of galaxies is usually assumed to be r2oo, indicating the radius for which 
the energy density of the cluster becomes 

P200 = 200p cr , where p cr = 4.6975 x 1CT 30 h 2 50 gr / cm 3 [32]. 

Using equation (72) we find 

r cr = 9lM^{^fh^M P c. (74) 
The total geometrical mass corresponding to this value is 

M{r) = 4.83 x 10 16 /3 3 / 2 {^f K^M Q , (75) 
which is consistent with the observational values for the virial mass of clusters [32] . 
5.2 Radial velocity dispersion in galactic clusters 

Radial velocity dispersion in galactic clusters plays an important role in estimating the virial mass of 
the clusters. It can be expressed in terms of the virial mass as [34] 

Mv = ^jR v . (76) 

Assuming that the velocity distribution in the cluster is isotropic, we have (v 2 ) = (v 2 ) + (vg) + {vfy = 
3(v 2 } = 3a 2 , the radial velocity dispersion a 2 for clusters in warped DGP model can be obtained 
from equation (37) as 

±{po*) + lp\' = 0, (77) 

where o\ and a r are related by 3of = a 2 . On the other hand, by neglecting the cosmological constant 
the Einstein field equation (43) becomes 

Integrating, we obtain 

where C\ is an integration constant. By eliminating A' from equations (77) and (79), we obtain 
f-, , X b 4 2\ d ( 2^ k\M(t) k\M(t) Ci 

I 1 + ^ ) Tr^ = -^^ P(r) " -^ 2 ~ f>{r) ~ V 2 ^ m 
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Figure 2: Left, The radial velocity dispersion in warped DGP brane-world model and right, the same parameter in 
DGP brane-world scenario for the NGC 5813 cluster with /3 = 0.766, r 20 o = 0.87Mpc, r c = 25Kpc, k B T g = 0.52A"eV 
and Ci = 0.503, 1.005 x 10 s , 2.011 x 10 8 M©, C 2 = 0.02,0.03,0.04 M%/Kpc 4 for solid, dashed and dot-dashed curves, 
respectively. 



Integration now gives the following solution 

_2 



where C2 is an integration constant. For most clusters /3 > | and therefore, in the limit r >> r c , the 
gas density profile (65) can be written as [23] 



, r \ ~ 3/3 2 
P,W = Po(-J > ( 82 ) 



Now, substituting equations (45), (73) and equation (82) into equation (81), for /3 / 1 we obtain 

1 + h 4 A * _ _ _ ^ v f^V w + *dk + ft i _ <± (L)» , (83) 



6 ^ J r 12(1-^(1-3/3) W /im p (l + 3/3)r Po 

and for /3 = 1 we find 

(l + *t«^) ^ = S^rl ( ± + ^)(L) 3 + *dk + ^ i - ^ f llV . (84) 
V 6 5H / ' 8 c Ur" r 4 / W /on, 4 r m \rj 

As is well-known, the simple form of(r) = B/(r + b) for the radial velocity dispersion and the relation 
p(r) = A/r(r + a) 2 for the density of the galaxies in cluster, with B, b, a and A constants, can be 
used to fit the observational data [34]. For r « a, p(r) ~ A/r, while for r » a, p(r) behaves like 
p(r) ~ A/r 3 . Here, our expression for of can be also used to fit the observational data. Therefore, 
the comparison of the observed velocity dispersion profiles of the galaxy clusters and the velocity 
dispersion profiles predicted by the warped DGP brane-world model may give a powerful method to 
discriminate between the different theoretical scenarios. 

Finally, we compare the radial velocity dispersion in warped DGP brane-world model with the 
radial velocity dispersion in the other theoretical models. As we noted before for A^ 5 ) = = 0, the 
warped DGP model reduce to the DGP model and equation (84) for (3 = 1 reduce to 

4 _ (i + ^) (z\> + Mi + ft I _ ft ( lY , , 85) 



4r 4 r 4 / \r c J pm p 4 r po 

which is the radial velocity dispersion in DGP model, equation (65) in [19]. The radial velocity 
dispersion of galaxy clusters in Palatini f(R) gravity for 7 = 3, which is corresponding to j3 = 1, is 
also presented as [18] 

r 2_ a f F ' -3, , k sT g ( 1 lnr\ cl d 3 
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where F(R) = d ^j^ . For f(R) = R this relation reduces to 

2 = ^/l h^X + ^ cW 

which is equation (84) with 8ttG = k\ and A^ 5 ** = A5 = 0. The same relation has been also obtained 
in Randall-Sundrum II model with this difference that the origin of geometrical mass in f(R) gravity 
is the extra terms in the Einstein-Hilbert action whereas in the latter is the global bulk effect. In 
figure 2 we have plotted the radial velocity dispersion for the cluster NGC 5813. The numerical 
values of it are in the ranges /3 = 0.766, r c = 25Kpc, ksT g = 0.52K ev, r2oo = 0.087Mpc [32] and the 
radial velocity is about 240km /s [35]. As one can see the radial velocity dispersion in warped DGP 
brane-world is compatible with the observed profiles and for the same value of constants C\ and C2 is 
slower than the DGP brane-world model. 



6 Conclusions 

The virial theorem plays an important role in astrophysics because of its generality and wide range of 
applications. One of the important results which can be obtained with the use of the virial theorem is 
to derive the mean density of astrophysical objects such as galaxy clusters and it can be used to predict 
the total mass of the clusters of galaxies. In the present paper, using the collision-less Boltzmann 
equation, we have obtained the generalized virial theorem within the context of the warped DGP 
brane-world model. The additional geometric terms due to the induced curvature term on the brane 
and non-local bulk effect in the modified gravitational field equations provide an effective contribution 
to the gravitational energy, equation (44) , which may be used to explain the well-known virial theorem 
mass discrepancy in clusters of galaxies. Finally, we have compared the virial theorem results with 
the observational data for galaxy cluster which can be obtained from the X-ray observation of the 
gas in the cluster and expressed the geometrical mass in term of observational quantities, like the 
temperature and the gas profile density. 
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